require(plotrix)
x<-rnorm(n=20,mean=7,sd=1.5)
x
y<-rnorm(n=20,mean=5,sd=1.0)
y
octopus<-function(x,y)
  
{
  X11()
  par(mfrow=c(3,2))
  multhist(list(x,y),col=c("green","red"),main="Peso de paralarvas \n de O.vulgaris aps ecloso",xlab="Peso (mg)",ylab="Dias aps ecloso")
  boxplot(list(x,y),main="Boxplot das amostras \n de paralarvas")
  qqnorm(x,xlab="Quantis esperados de x",ylab="Quantis das amostras x",main="Quantil esperado de x e \n Quantil das amostras de x")
  qqline(x)
  qqnorm(y,xlab="Quantis esperados de y",ylab="Quantis das amostras de y",main="Quantil esperado de y e \n Quantil das amostras de y")
  qqline(y)
  plot(y~x,xlab="Amostra x",ylab="Amostra y", main="Amostras x e y")
  par(mfrow=c(1,1))
  octopus.amostra1<-summary(x)
  octopus.amostra2<-summary(y)
  octopus.coef<-coef(lm(y~x))
  return(list(octopus.amostra1,octopus.amostra2,octopus.coef))
}

#teste1
octopus(x,y) # DEU CERTO!!!\o/
#"-Crebro, o que vc quer fazer esta noite?" 
#"-A mesma coisa que fazemos todas as noites: tentar conquistar o mundo!!!"
